Proton-induced magnetic order in carbon: SQUID measurements 1 

J. Barzola-Quiquia, S. Petriconi, P. Esquinazi*, M. Rothermel, D. Spemann, A. Setzer, 

T. Butz 

Institut fur Experimentelle Physik II, Universitdt Leipzig, Linnestrasse 5, D- 04 103 Leipzig, Germany 



Abstract 

In this work we have studied systematically the changes in the magnetic behavior of highly oriented pyrolytic graphite 
(HOPG) samples after proton irradiation in the MeV energy range. Superconducting quantum interferometer device 
(SQUID) results obtained from samples with thousands of localized spots of micrometer size as well on samples 
irradiated with a broad beam confirm previously reported results. Both, the para- and ferromagnetic contributions 
depend strongly on the irradiation details. The results indicate that the magnetic moment at saturation of spots of 
micrometer size is of the order of 10 -10 emu. 
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1. Introduction 

Magnetic order at room temperature in metal-free 
carbon-based structures remains one of the exciting 
issues in fundamental and applied research across 
all scientific disciplines. However, the lack of repro- 
ducibility of early results added to the unknown, 
in some cases late characterization of the magnetic 
impurities [1] increased substantially the scepticism 
of the scientific community. In the year 2003 some 
of us reported that proton irradiation of MeV en- 
ergy on HOPG samples triggers ferro- or ferrimag- 
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netism at room temperature [2,3]. Although reports 
on the reproducibility of this behavior in two differ- 
ent HOPG samples followed the original publication 
[4] , the lack of independently published studies that 
successfully triggered magnetic order after ion irra- 
diation of carbon pushed us to elaborate and extend 
our experimental work. 

Although there are several theoretical works in 
the literature on magnetic order in carbon, specially 
hydrogen-induced (see for example Ref. [5] and 
references therein), apart from the proton induced 
magnetic order in graphite there are no much new 
systematic experimental works that show ferromag- 
netism in metal-free carbon, including in particular 
a rigorous characterization of the samples impuri- 
ties. The study done in Ref. [6] reported that Nitro- 
gen and Carbon irradiation of nanosized diamond 
powder triggers magnetic order at room tempera- 
ture. However, no study of the impurity concentra- 
tion was presented in that work. The studies done 
in [7] reveal that carbon films prepared by CVD 
on stainless steel substrates reach magnetization 
values of the order of 0.15 emu/g at room tem- 
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perature, comparable to those reported in earlier 
studies [8,9]. In that work the amount of measured 
impurities appears to be not enough to account for 
the absolute value of the magnetic moments of the 
samples. In Ref. [10] the magnetization of proton 
irradiated graphite was studied after one low-dose 
irradiation. The authors show there that proton 
irradiation induces Curie-type paramagnetism in 
graphite but no measurements of the hysteresis 
loops before and after irradiation were apparently 
done in order to check for the induced magnetic or- 
der. The experimental difficulties to reproduce this 
magnetic order and the weakness of the ferromag- 
netic signals, which are sometimes at the limit of 
the sensitivity of current experimental characteri- 
zation methods, including SQUID's, are the main 
obstacles that preclude a rush development of this 
interesting and important subject. 

The aim of this report is threefold. Firstly, we dis- 
cuss the expected ferromagnetic signals after pro- 
ton irradiation based on our previous publication 
[2] and the necessary precautions one needs to take 
for the SQUID measurements of these weak signals. 
To increase the relative sensitivity of these measure- 
ments and to avoid experimental artifacts as well 
as any introduction of magnetic impurities through 
sample handling, we developed a new, simple sam- 
ple holder that allows consecutive irradiation and 
SQUID measurements without touching the sam- 
ple or its holder as well as their relative positions. 
Secondly, in order to check the magnetic impurity 
concentration of the HOPG samples and the accu- 
racy of the method we used (Particle Induced X- 
ray Emission (PIXE)), we have compared the impu- 
rity concentrations measured with three experimen- 
tal methods on similar HOPG samples. Thirdly, we 
have irradiated three HOPG samples under different 
conditions and measured the changes in their mag- 
netic behavior with a SQUID. We show that proton 
irradiation induces ferromagnetism in graphite, as 
demonstrated earlier [2], and we extend our studies 
to the induced paramagnetism. Through the irradi- 
ation of thousands of spots of micrometer size in a 
single sample we were able to measure their ferro- 
magnetic signal and estimate the average magnetic 
moment of each spot. 

2. Experimental details 

The studies performed in Ref. [2] indicate that 
the ferromagnetic magnetic moment at saturation 



m s depends on the total implanted charge C t . 
Roughly speaking m s ~ 2 x 10~ 7 [emu/^C°- 5 ]C t - 5 
for C t < 10 3 /iC (the usual magnetic moment (cgs) 
unit "emu" is equal to 10~ 3 Am 2 ). C t means the to- 
tal charge irradiated on the same area. The square 
root dependence suggests that at small doses the 
experimental signal should be directly proportional 
to the irradiation dose, which produces vacancies, 
adatom defects and/or vacancy- hydrogen com- 
plexes. However, it tends to saturate at higher 
doses, probably due to damage accumulation (e.g. 
saturation of dangling bonds), hydrogen outgassing 
and/or annealing of certain defects [11]. This square 
root relation should be taken only as a rough es- 
timate for the expected ferromagnetic signal and 
for broad irradiation areas. The actual value of m s , 
however, depends on parameters like proton cur- 
rent, i.e. the higher the current the lower is m s , 
(e.g. due to heating effects), or beam size (broad or 
narrow beam irradiation) [4]. According to earlier 
results [2] we expect m s of the order of 10~ 6 emu 
for a total irradiated charge of 100 /iC. That means 
that the SQUID magnetometer should provide re- 
liable and reproducible results within a magnetic 
moment range better than 1 /xemu, specially after 
introducing and taking out the sample with holder 
several times into or from the SQUID. 

Usually at fields B < 1 T the reproducibility 
of commercial SQUID's is better than 1 /xemu and 
therefore these magnetometers can be used to mea- 
sure the effects produced by irradiation (an exam- 
ple is shown below). Care should be taken, how- 
ever, with possible artifacts of these systems, spe- 
cially magnetic field hysteresis due to the electron- 
ics and/or superconducting solenoid properties [12]. 
The reproducibility of each SQUID system should 
be checked before starting the irradiation steps. 

The magnetic moment measurements were per- 
formed with a SQUID magnetometer from Quan- 
tum Design with the reciprocating sample option 
(RSO). We note that the SQUID sensitivity without 
this option is not enough to measure accurately the 
effects produced by irradiation, specially when the 
magnetic signal is of the order of /xemu for the sat- 
uration ferromagnetic moment. The magnetic field 
was applied parallel to the graphene planes in all 
measurements in order to diminish the diamagnetic 
background. As an example for the reproducibil- 
ity and error, Fig. 1 shows the difference between 
hysteresis loops of the same HOPG sample (includ- 
ing holder) before irradiation measured at different 
days. Each measurement was performed following 
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Fig. 1. Difference of magnetic moments as a function of field 
(hysteresis loops) measured for the same HOPG sample and 
holder (in this case number 3 as example) at four different 
days at T = 100 K. The difference Am is calculated tak- 
ing the measurement at a certain day as reference. Similar 
differences are obtained at all temperatures and choosing 
measurements at other days as reference. 

the same sequence after removing and introducing 
the whole sample with holder from/into the SQUID 
apparatus, leaving the sample-holder several days 
at ambient conditions between the measurements. 
The difference Am plotted in Fig. 1 should be ide- 
ally zero at all fields. We recognize, however, that 
both paramagnetic- as diamagnetic-like deviations 
are obtained after subtraction of an arbitrarily cho- 
sen hysteresis loop (measured at a certain day) from 
the other loops. The largest deviation obtained is of 
the order of 8 x 10 -7 emu at a field of 1 T. Because 
most of the HOPG samples we have measured in 
their virgin states show ferromagnetic-like hystere- 
sis loops, which origin is not related to the magnetic 
impurities [13], the effects produced by the irradia- 
tion are much better identified doing this point-by- 
point (at equal (B, T)) difference in Am. 

Because we do not know the amount of sample 
that remains magnetic after irradiation, all the 
SQUID data are presented as magnetic moment 
directly measured with this method. Although one 
knows the penetration depth of protons of the used 
energy in graphite (< 48 /xm) it is still unclear, 
which is the range where the main ferromagnetic 
and paramagnetic signals come from. Therefore, 
it has little sense to divide the measured values of 
magnetic moment by the total sample mass. Taking 
into account recently done x-ray magnetic circular 
dichroism and magnetic force microscopy measure- 
ments of irradiated spots in 200 nm carbon films, 



which provide clear evidence for the existence of 
magnetic order at the spot position [14], we tend to 
assume that the ferromagnetic layer in our samples 
should be mainly located at the first micrometer 
from the sample surface. 

The irradiations were done with the high-energy 
nanoprobc LIPSION, a single ended 3MV SIN- 
GLETRON accelerator with a RF-source for pro- 
tons and alpha particles. Two irradiation chambers 
allow the irradiation with proton beams of diam- 
eters as low as 50 nm up to 0.8 mm at different 
proton currents at MeV energy 

All three samples were HOPG from Advanced Ce- 
ramics with a rocking curve width of 0.4° (grade A), 
-(a) Sample 1 had a mass of 3.5 mg with a size 
2 x 3 x 0.3 mm 3 . As in Ref. [2], it was glued on a 
high purity Si substrate with varnish. The irradi- 
ation of this sample consisted on three broad pro- 
ton spots of 2 MeV proton energy, 0.8 mm diameter 
each, 0.7 mm distance, 150 /iC total charge per spot 
and a proton current of 100 nA. 
-(b) Sample 2 of 5.7 mg weight and 2 x 4 x 0.3 mm 3 
size was glued in the middle of a 13.5 cm long and 
1 mm diameter pure Cu wire. Due to the long length 
and the homogeneous magnetic response of the wire, 
it provides a small contribution to the SQUID sig- 
nal. The whole ensemble sample plus wire including 
the corresponding couplings for the SQUID and ac- 
celerator chamber, were used for the irradiation of: 
first 10 4 spots and later 2 x 10 4 spots more. Each 
spot was made with a 2.25 MeV proton energy, beam 
diameter of ~ 2 /um, a total charge of 1.16 nC (flu- 
ence 0.37 nC/^m 2 = 2.3 x 10 17 protons/cm 2 ) and 
we used 700 pA proton current. 
-(c) Sample 3 had a mass of 10.8 mg and a size of 
4 x 3 x 0.4 mm 3 . It was glued with varnish on the 
middle of a 13.5 cm long pure quartz rod. Before fix- 
ing it a 20 nm gold film was deposited on the rod, 
which is necessary for the measurement of the proton 
current in the accelerator chamber during irradia- 
tion. For the irradiation of spots we used a ~ 1.5 /zm 
diameter proton beam, which irradiated a ~ 5 yum 
diameter area following a computer controlled spi- 
ral movement. 400 spots (25 /im apart) with a total 
charge of 191 /uC, fluence per spot of 24.3 nC/^m 2 
(1.52 x 10 19 cm' 2 ) at a proton current of 6.5 nA 
(10 times larger than for sample 2) and 2.25 MeV 
energy were prepared in sample 3. 
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3. Results 



3.1. Impurity measurements 

As described in section 1, we expect ferromagnetic- 
like signals of the order of a few /xemu. How much 
ferromagnetic Fe is necessary to produce a magnetic 
moment of, e.g., 5 /jemu? From literature data one 
estimates easily that 23 ng pure Fe (or a volume of 
~ 3 x 10~ 9 cm 3 ) would be enough to produce this 
magnetic moment assuming that this small amount 
is ferromagnetic at room temperature. Under these 
assumptions the relative Fe concentration in a typ- 
ical HOPG sample would be of the order of 6 ng/g- 
That means that we need impurity measurements 
that provide a sensitivity of at least 1 ng/g for Fe. 
The method called Particle Induced X-ray Emis- 
sion (PIXE) has this sensitivity (or better) for the 
analysis of Fe in a carbon matrix. Our PIXE mea- 
surements performed in situ and during irradiation 
show that the total amount of Fe impurities in our 
HOPG samples is 0.6 ± 0.04 /ig/g (i.e. 0.15 ppm)[4] 
for grade A samples. To check the accuracy of our 
PIXE analysis we have measured similar HOPG 
samples with two other methods. Neutron Activa- 
tion Analysis (NAA) on a 30 mg HOPG grade B 
sample with an neutron activation time of 18 days 
and 7-ray measuring time of 6 hours provides a total 
Fe impurities of 0.17 ± 0.03 /ig/g. Our PIXE mea- 
surements on the same sample gave for Fe 0.17 (J-g/g 
as well. The third method used was x-ray fluores- 
cence with a EDXRF-spectrometer Quan X. The Fe 
concentration in the same HOPG grade A sample 
measured with PIXE was below the minimum de- 
tection limit (~ 5 /ig/g) this analytical method has. 
The amount of other magnetic impurities was much 
below that for Fe (for example, 2.5 x 10 -3 /xg/g for 
Co). 

Concluding, assuming the worst, unlikely case, 
the maximum ferromagnetic magnetic moment at 
saturation one expects from this amount of mag- 
netic impurities in our HOPG samples - were this 
amount ferromagnetic at room temperature - is < 
5 x 10~ 7 emu. Certainly, one should rule out that 
due to an improper sample handling and between ir- 
radiation steps a small Fe grain with a mass of a few 
tens of ng does get fixed somewhere at the surface 
of the HOPG sample. Therefore, special holders as 
well as systematic irradiation steps are necessary. 
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Fig. 2. (a) Total magnetic moment (HOPG sample 1 with 
the Si substrate) as a function of temperature at a constant 
magnetic field applied parallel to the graphene planes for the 
sample before and after proton irradiation with a broad beam 
(total charge 450 fiC at 100 nA proton current), (b) The 
difference between the two curves from (a). This difference 
reveals directly the irradiation effect. The continuous line is 
the function 3 X 10 _5 [emu K]/T + 3 X 10 _6 [emu]. 

3.2. Broad irradiation 

The temperature dependence of the magnetic mo- 
ment before and after irradiation of sample 1 is 
shown in Fig. 2(a) in a semilogarithmic scale. The 
temperature dependence of the virgin curve shows a 
minimum (maximum diamagnetism) at T ~ 30 K, 
which is usual for HOPG samples of good qual- 
ity. After a broad proton beam irradiation cover- 
ing most of the sample area, the magnetic moment 
shows a clear increase in all the temperature range. 
Figure 2(b) shows the difference between the mag- 
netic moment after the irradiation minus that of the 
virgin state as a function of temperature at a con- 
stant field of 1 T. This difference can be roughly un- 
derstood as the sum of two contributions, namely, a 
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paramagnetic one, which follows roughly the Curie 
law 3 x 10~ 5 /T emu, and a ferromagnetic constant 
contribution 3 x 1CT 6 emu, i.e. m(T,B > B s ) ~ 
h 3 x 10~ 6 , where B s is the minimum satura- 
tion field for the ferromagnetic part. The small but 
clear deviation between the fit and the data shown 
in Fig. 2(b) can be interpreted as follows. The main 
variable of the Brillouin function Bj(x) with x = 
gJ/iBB/ksT reaches relatively high values at low 
temperatures for an applied field of 1 T. Assuming 
for simplicity the product of the Lande factor g with 
the total angular momentum J, gj ~ 1, this vari- 
able is 0.13 < x < 2.2 x 1(T 3 for 5 K< T < 300 K 
at B = 1 T. Only for x <C 1 one is allowed to keep 
the first term of the Brillouin function that provides 
the simple 1/T Curie law. Part of the deviation may 
also come from the assumption of a strictly temper- 
ature independent ferromagnetic contribution. 

The hysteresis loops shown in Fig. 3, obtained at 
two temperatures subtracting the loops after irradi- 
ation from those obtained in the virgin state, justify 
the assumption of the two magnetic contributions. 
The increase of the magnetic moment after irradia- 
tion at room temperature is mainly due to the in- 
crease of the ferromagnetism of the sample; only <~ 
25% of the increase is due to a paramagnetic contri- 
bution at 300 K and 1 T. The paramagnetic contri- 
bution is clearly recognized in Fig. 3 from the slope 
of the loops at fields above <~ 0.25 T. The inset in 
this figure shows clearly the finite irreversibility pro- 
duced by the irradiation with coercivity fields of the 
order of 0.02 T. 

Taking into account the total irradiated charge 
per spot and after Ref. [2] we expect a m s (300) <~ 
3 x 2 x 10~ 7 Vl50 ~ 7 x 10 -6 emu, in comparison we 
obtain m s (300) ~ 3.0xl0~ 6 emu. The results shown 
in Fig. 2 clearly indicate that broad irradiation - at 
the used proton current and fluence - triggers two 
magnetic contributions, one due to independent, lo- 
calized magnetic moments (e.g. dangling bonds due 
to the disorder produced by the irradiation) and a 
second one with all the characteristics of magnetic 
order with a Curie temperature above room tem- 
perature. We stress that the effect on the magnetic 
properties of graphite due to proton irradiation de- 
pends on several parameters as the total implanted 
or irradiated charge, fluence and proton current as 
well as on the geometry of the used proton beam, as 
the next section discusses. 
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Fig. 3. Hysteresis loops for sample 1 obtained from the dif- 
ference between the loops measured after and before irradi- 
ation Am = m a — rrtj, at the same magnetic fields in each 
state at two temperatures. The inset blows up the data in a 
smaller field range. 

3.3. Magnetic spots of micrometer size 

3.3.1. Low proton current 

Localized proton irradiation of spots of micro- 
meter size triggers a slightly different magnetic re- 
sponse on graphite. Figure 4 shows the total mag- 
netic moment (sample and holder) as a function of 
temperature at a field of 1 T for sample 2 in three 
different states, virgin, with 10 4 and with 3 x 10 4 ir- 
radiated spots, each of them with a total charge of 
1.16 nC. The temperature dependence of m is quali- 
tatively similar to that of Fig. 2. The difference is 
mainly due to the sample misalignment respect to 
the applied field and partially also to a different con- 
tribution of the sample holder. 

The irradiation effect on the temperature depen- 
dence can be better recognized from the difference 
between the m(T) curves shown in the inset of 
Fig. 4. The inset shows that no Curie-like (1/T)- 
contribution is obtained after producing the local- 
ized spots. Part of the magnetic response of the 
spots is due to ferromagnetism, as the hysteresis 
loops shown in Fig. 5 indicate. From these loops 
we obtain a saturation magnetic moment m s ~ 
10~ 6 emu. For the first irradiation with only 10 4 
spots the error involved in the subtraction, specially 
at low fields, does not permit to assure the existence 
of an hysteresis. We stress that signals of the order 
of 2 x 10~ 7 emu are at the limit of reliability (not 
resolution!) of our SQUID apparatus. 

It is well known that the paramagnetic Curie law 
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Fig. 4. Total magnetic moment of sample 2 as a function of 
temperature at a field of 1 T for the virgin (•), after the first 
irradiation with 10 4 spots (total charge 11.6 /iC at 700 pA 
(□)), and after the second irradiation adding 2 X 10 4 similar 
spots as in the first irradiation (if). The inset shows the 
difference of magnetic moment between the first (I) (second 
(*)) irradiation and the virgin state. The continuous lines are 
fits to the function mo +rai exp(— T/Tq) with the parameters 
m = 2.2(5.6) X 10~ 6 emu, mi = 1.9(5.07) X 10" 6 emu and 
To = 116(140) K for the first (second) irradiated sample. 

holds only if x <C 1. This law is a consequence of 
thermal average involving (2 J + 1) equally spaced 
levels, which originate from the effect of the applied 
field on one multiplet. If at large magnetic fields new 
multiplet levels start to contribute to the statisti- 
cal average or if their energy levels are not equally 
spaced, then deviations of m(T) from the Curie law 
to a weaker T— dependence are expected. The ob- 
served positive curvature for fields above 0.25 T in 
the hysteretic loops, see inset in Fig. 5, suggests that 
the applied field influences the number of multiplet 
levels that determine the total magnetic moment. 
The clearly weak T— dependence of the irradiated 
sample 2, see inset in Fig. 4, agrees with this expecta- 
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Fig. 5. Hysteresis loops for sample 2 obtained from the dif- 
ference between the loops measured after and before irra- 
diation Am = m a — mi, at the same magnetic fields in the 
states after two irradiations at 100 K. The points (o) were 
obtained after irradiation of 10 4 spots (11.6 /iC total charge 
at 700 pA proton current) and (•) adding 2 X 10 4 (34.8 fiC 
total charge) spots on a different area of the same sample. 
The bars indicate the maximum expected error due to the 
reproducibility of our SQUID and the subtraction. The in- 
set shows the loops in a broader field range after the first 
(dashed line) and second irradiation (continuous line). 

tion. The T— dependence shown in this inset can be 
fitted with an exponential function of the form m ~ 
mo + mi exp(— T/Tq) with mo,i and To free fitting 
parameters. If we use the Brillouin function plus a 
T— independent ferromagnetic contribution we can 
in principle fit the measured T— dependence but us- 
ing gj values that are above any reasonable limit, 
indicating its inadequacy to understand the mag- 
netism of the irradiated spots. Whatever the rea- 
sons for the observed behavior, the measured curves 
indicate a constant temperature ferromagnetic-like 
term of the order of mo ~ 2x 10 -6 emu for the sam- 
ple with the first 10 4 spots and mo — 5.6 x 10~ 6 emu 
for the sample with 3 x 10 4 spots at 1 T. These 
numbers indicate a magnetic moment of the order of 
2 x 10~ 10 emu per spot produced at the conditions 
described in section 2. 



3.3.2. Large proton current 

From MFM (Magnetic Force Microscopy) mea- 
surements we know that higher currents decrease the 
magnetic phase contrast at the center of the spot, 
indicating the vanishing of magnetic order in part of 
the irradiated area [4] . To test this behavior with the 
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Fig. 6. Temperature dependence of the difference between the 
magnetic moments measured after and before the irradiation 
of 400 micrometer-size spots with a total charge of 191 fiC 
at a current of 6.5 nA in sample 3. The measurements were 
done at a constant field of 0.3 T. The continuous line follows 
the equation 3.9 x 10~ 7 [emu] +1.84 x 10~ 6 [emu K]/T. 



SQUID the spots in sample 3 were produced with a 
<~ 9 times larger proton current than for the spots 
produced in sample 2. Figure 6, as Fig. 2(b), shows 
the temperature dependence of the magnetic mo- 
ment produced by the irradiation of 400 spots (flu- 
ence 24.3 nC/^m 2 = 1.52 x 10 19 cm~ 2 , total charge 
191 /j,C). In contrast to the effects of the spots pro- 
duced in sample 2, the irradiation in sample 3 trig- 
gers a paramagnetic contribution that follows the 
Curie law plus a smaller, constant ferromagnetic- 
like background. This small, last contribution is rec- 
ognized also in the s-shape curve at low fields ob- 
tained from the difference between the hysteresis 
loops. The results indicate that it is not only the to- 
tal implanted charge what determines the ferromag- 
netic and/or the paramagnetic effect but also the 
proton current appears to have an important role 
due to overheating effects. 



We have performed high sensitive measurements 
with a commercial SQUID with RSO option of the 
effects produced by proton irradiation in the MeV 
energy range on HOPG samples. All proton irradi- 
ations produce two magnetic contributions to the 
magnetic properties of HOPG. One contribution de- 
pends clearly on temperature and in some cases fol- 
lows a Curie-like behavior. The second one can be 
attributed to ferro- or ferrimagnctism. These con- 
tributions are supported by the measured hysteresis 
and s-form of the loops as a function of magnetic 
field. The details of these contributions depend not 
only on the total implanted charge, the fluence but 
also on the proton current and the beam geome- 
try, i.e. broad beam irradiation or localized spots. 
Largest ferromagnetic signals were obtained for lo- 
calized spots produced at relatively low proton cur- 
rents (< 1 nA) and ffuences < 1 nC//im 2 (6.25 x 
10 cm~ 2 ) . The saturation magnetic moment per 
micrometer spot is of the order of 10 -10 emu. The 
absolute values of the observed effects are of the 
same order as those published in the original work 
of Ref. [2] . The observed decrease of the effects at 
large fluences and proton currents agrees with the 
decrease of the (magnetic) phase contrast with flu- 
ence and current observed in MFM measurements 
on similar spots produced in HOPG surfaces [4]. 

We note that both, SQUID and MFM measure- 
ments of the magnetic behavior of micrometer size 
spots indicate that upon irradiation conditions they 
may not behave ferromagnetically with a finite re- 
manent magnetic moment at zero field. In this case 
the characterization of the spots by means of MFM 
under zero field conditions is difficult since no sig- 
nificant contrast difference will be measured after 
a change of the tip polarization or after applying a 
magnetic field on the spots. 

Three experimental methods were used to check 
the amount of magnetic impurities in similar HOPG 
samples from the same company. The results indi- 
cate that the total amount of magnetic impurities 
is much below that needed to understand quanti- 
tatively the measured ferromagnetism after irradi- 
ation. With the implementation of a special sam- 
ple holder for SQUID measurements and irradiation 
runs, we could increase the reliability of the SQUID 
results and rule out the influence of artifacts due to 
handling on the sample and its holder. 

On the origin of the magnetic order in proton irra- 
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diated carbon we note the following. The overall re- 
sults suggest that carbon defects (e.g. adatoms or va- 
cancies) should play a role in the magnetic order ob- 
served. According to recent models [11] hydrogen at 
carbon defects (H- vacancy and -adatom complexes) 
may also contribute triggering or even enhancing a 
local magnetic moment. Regarding the role of hy- 
drogen and taking into account that HOPG samples 
have a substantial amount of hydrogen at the first 
micrometer from the surface [15], one can speculate 
that not the implanted charge but the dissociation 
of molecular hydrogen (already in the sample) pro- 
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